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Abstract. We test the present expansion of the universe using supernova type la 
data without making any assumptions about the matter and energy content of the 
universe or about the parameterization of the deceleration parameter. We assume the 
cosmological principle to apply in a strict sense. The result strongly depends on the 
data set, the light-curve fitting method and the calibration of the absolute magnitude 
used for the test, indicating strong systematic errors. Nevertheless, in a spatially flat 
universe there is at least a 5a evidence for acceleration which drops to l.Scr in an open 
universe. 
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1. Introduction 

The first evidence for an accelerated expansion of the universe at the present epoch 
was presented in the late 1990s by Riess et al. [1] and Perlmutter et al. [2]. Fitting a 
cosmological model including a cosmological constant matched the supernovae type la 
of the considered data sets significantly better than a model without such a constant. 
In the following years, a lot of effort was put into getting larger supernova (SN) sets of 
increasing quality [3, 4, 5, 6]. On the other hand, a variety of cosmological models have 
been developed that could be fit to the newest data sets. Those models are characterized 
by certain parameters which are used to derive a redshift-luminosity relation that can be 
compared to the observed values. The problem with this so-called dynamical approach 
is that it is impossible to test without assumptions on the matter and energy content 
of the universe if there really is a phase of acceleration in the expansion history of the 
universe. 

Some authors tried to avoid this problem by taking a kinematical approach, i.e. 
they only considered the scale factor a and its derivatives, such as the deceleration 
parameter q{z), without using any model specific density parameters or a dark energy 
equation of state. The first kinematical analysis of SN data was done by Turner and 
Riess [7] who considered averaged values qi for redshift z < zi and q2 for z > zi 
concluding that a present acceleration and a past deceleration is favoured by the data. 
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Other authors tested a variety of special parameterizations of q{z) [8, 9] or used a{t) 
[10] or the Hubble rate H{z) [11]. Instead of considering a special parameterization, a 
more general approach has been done by Shapiro and Turner [12] who expanded the 
deceleration parameter q into principal components. Rapetti et al. [13] expanded the 
jerk parameter j into a series of orthonormal functions. However, one has to be careful 
when doing a series expansion in z as SNe with a redshift larger than 1 are not within 
the radius of convergence. This problem can be solved by reparameterizing redshift [14]. 

In the present work we will neither make any assumptions about the content of 
the universe, nor about the parameterization of q{z) or another kinematical quantity. 
Moreover, we do not need to assume the validity of Einstein's equations. Instead, we 
will only ask the question if the hypothesis holds that the universe never expanded 
accelerated between the time when the light was emitted from a SN and today. Our 
assumption for the test is that the universe is isotropic and homogeneous. This is 
certainly not true for the real universe in a strict sense. However, we follow the 
standard approach in assuming that the cosmic structure does not modify the observed 
SN magnitudes and redshifts apart from random peculiar motion. The basic idea of 
this analysis has already been presented by Visser [15] and has been applied to SN 
data by different groups [16, 17]. But while Santos et al. [16] made some mistakes in 
their analysis (that we will discuss later in section 6) and Gong et al. [17] only state 
that accelerated expansion is "evident" , we are able to give a quantitative value for this 
evidence. Additionally, we study the size of systematic effects. 

For the fit one usually marginalizes over a function of the absolute magnitude M and 
the Hubble constant Hq because these two values cannot be determined independently 
by only considering SNe. Marginalization is not suitable for our analysis because in 
order to do so a special cosmological model or at least a parametization of q{z) has to 
be inserted, which is what we want to avoid. But the SNe can be calibrated by cepheid 
meassurements and thus the values of M and Hq are determined. Using this additional 
information, no marginalization is needed. As there is still a controversy [18] about the 
appropriate calibration method, we will take two quite different results of calibration 
into account [19, 20]. 

In order to test the robustness of our analysis, we consider two different SN data sets 
(the 2007 Gold sample [5] and the ESSENCE set [6]), where the data of the ESSENCE 
set are once obtained by the multicolour hght-curve shape (MLCS2k2) [3] fitting method 
and once by the spectral adaptive hght-curve template (SALT) [21]. Both calibration 
methods are applied to each set. The analysis shows that in all cases the data indicate 
an accelerated expansion. But the confidence level at which acceleration can be stated 
strongly depends on the data set, the fitting method and the calibration. We begin our 
analysis by assuming a flat universe, but later also consider the cases of an open and a 
closed universe. 
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2. Method 

We want to keep our test as model-independent as possible, but still a few assumptions 
have to be made. We consider inflationary cosmology to be correct. This implies large 
scale homogeneity and isotropy of the universe as well as spatial flatness. (Yet, we 
will give up the assumption of a flat universe later in section 5.) We also assume that 
cosmological averaging (here along the line of sight) does not modify the result obtained 
in a Friedmann-Lemaitre model. In such a universe, the luminosity distance di,{z) is 
given by 

As we are interested in the question whether the universe accelerates or decelerates at 
the present epoch, we need to examine the deceleration parameter q{z): If q is positive, 
the universe expands decelerated, for a negative q it accelerates. q{z) can be expressed 
in terms of the Hubble parameter H{z): 

q{z)^^{l + z)-l, (2) 

where the prime denotes the derivative with respect to z. Integrating this equation 
yields 

Ho Jo l + z 

Our null hypothesis is that the universe never expanded accelerated, i.e. q{z) > 
for all z. Under this assumption, the above equation turns into the inequality 
H(z) /"^ 1 

ln-^> / -—dz^Hi + z) (4) 

or H{z) > Ho{l + z). Thus, for the luminosity distance we have 

1 /"^ d5 1 
dj.(z)<(l + z)— ^(l + z)—Hl + z). (5) 

-no Jo + ^ -"0 

In order to test our hypothesis, we consider different data sets of SNe type la. If 
the observed luminosity distance is signiflcantly larger than the luminosity distance of a 
universe with a constant q — 0, the hypothesis can be rejected at a high confldence level. 
Note that the rejection of the hypothesis does not mean that there was no deceleration 
between the time the light was emitted from a SN until now, it only gives evidence that 
there was at some time a phase of acceleration. Thus, we cannot determine the transition 
redshift between deceleration and acceleration. This restriction to our analysis comes 
from the integral over redshift in the calculation of di^{z). 

As a first step, the data in the considered sets have to be calibrated consistently. 
The distance modulus /j, is related to the luminosity distance by 

^ = m- M = 51ogrfL + 25, (6) 

where c^l is given in units of Mpc. The distance modulus of the SNe la in those sets is 
always given in arbitrary units because the absolute magnitude M cannot be measured 
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independently of the Hubble constant Hq. Only the apparent magnitude m and thus 
the relative distance moduli are measured at high precision. In order to determine the 
absolute magnitude, the SNe have to be calibrated by measuring the distance of cepheids 
in the host galaxies. Then also the Hubble constant can be determined. But there is still 
disagreement between different groups about the correct calibration analysis. In this 
work we will consider two results for M and Hq, namely that of Riess et al. [19] and that 
of Sandage et al. [20] . In the following we will refer to those results as Riess calibration 
and Sandage calibration, respectively. The difference in the calibration comes mainly 
from the different assumptions on the evolution of the cepheid period-luminosity relation 
with host metallicity [18]. So for preparing the data for our analysis, first the distance 
moduli fii have to be adjusted to the assumed absolute magnitude. 

We define the magnitude Afii as the observed distance modulus /ij of the i-th SN 
minus the distance modulus in a universe with constant deceleration parameter q = 
at the same redshift zf 

A/ii ^ - ij,{q ^ 0) ^ Hi - 5 log 



,1 + Zi) ln(l + Zi) 

-no 



- 25 . (7) 



If the error in redshift and the peculiar velocities of the SNe are already included in the 
error cr^. given in a certain data set, then the error cTj of A/Xj equals cr^.. Otherwise, the 
resulting error of A/ii is calculated by: 



2.1. ln(l + Zi) + 1 \ / 2 , 2\ 



Mi 



(8) 



Let Hth{z) be the theoretical distance modulus of a cosmological model, which 
describes the expansion of the universe correctly. Then 

Aiithizi) = A*th(^i) - = 0) (9) 

would be the "true" value corresponding to the meassured SN value A//j. The null 
hypothesis for our test is that the universe never expanded accelerated, i.e. A//th < for 
each SN. We reject that hypothesis if the measured value A/ij lies above a certain action 
limit y4a, otherwise we accept it. We want to keep the risk low that we conclude a late 
time acceleration of the universe, when there is indeed no acceleration at all. Therefore 
the action limit must be relatively high. If we want the confidence level for concluding 
an accelerated expansion to be 99%, the action limit must be Aa(99%) = 2.326cri. For 
a confidence level of 95%, the limit is As,{99%) ~ 1.645(7j. Here we assume that the 
measured values A//j at a given redshift follow a normal distribution. 

On the other hand, we can test the hypothesis that the universe expanded 
accelerated all the time from light emission of a SN until today. This hypothesis can be 
rejected at a 99% CL if A/Xj is below the action limit Ad{99%) = ^2.326(jj. That would 
mean that there must have been a phase of deceleration in the late time universe, but 
it would not exclude a phase of acceleration. 
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3. Data Sets 

The important parameters of the SNe are obtained by fitting their hght-curves. The 
results depend on the fitter that is used. The most common fitter is MLCS2k2 [3]. Here 
we will also consider the SALT fitting method [21]. The main difference between the 
two methods is in how the peak luminosity corrections are determined by the colour 
of the SNe [22]. While MLCS2k2 assumes that the corrections that have to be made 
are only due to dust, SALT takes an empirically approach to determine the relation 
between the colour and the luminosity correction. 

For the analysis we take the following SN la data sets: 

Gold 2007 (MLCS2k2): the Gold sample by Riess et al. (2007) [5], which was obtained 
by using the MLCS2k2 fitting method 

ESSENCE (MLCS2k2): the set given by Wood-Vasey et al. (2007) [6], which includes 
data from ESSENGE, SNLS [4] and nearby SNe [3], fitted with MLGS2k2 

ESSENGE (SALT): the same set fitted with SALT 

As suggested by Riess et al. [5], we discarded all SNe with a redshift oi z < 0.0233 
from the Gold sample, which leaves us 182 SNe with 0.0233 < z < 1.755. In ESSENGE 
(MLGS2k2) and ESSENGE (SALT), the SNe with bad fight curve fits were rejected for 
each fitting method scperately. This leaves 162 SNe for ESSENGE (MLGS2k2) and 178 
for ESSENGE (SALT) with 0.015 < z < 1.01. 153 of the SNe arc contained in both sets. 
Due to the differences of the SNe contained in ESSENGE (MLGS2k2) and ESSENGE 
(SALT) we will refer to them as different sets in the following. Nevertheless, you should 
keep in mind that they share a large number of SNe and thus are not independent sets. 

The Riess calibration yields a value for the V-band magnitude of My{to) — 
—19.17 ± 0.07mag, where to is the time of the B-band maximum. For the Gold 2007 
set Mv(to) "was considered to be — 19.44mag. Therefore, in order to get the appropriate 
distance modulus ^ = m — M, we need to substract 0.27mag from the value given in 
the Gold sample. From the distance modulus values in the ESSENGE sets 0.22mag 
have to be substractcd because in these sets a B-band magnitude of Mb = — 19.5mag is 
assumed and Riess et al. [19] give the relation Mb — My = —0.11. For this calibration 
one gets a Hubble constant of i^o = 73 ± 4 (statistical) ± 5 (systematic). 

The results of Sandage et al. [20] for the SNe absolute magnitudes are My — 
— 19.46mag and Mb = — 19.49mag and for the Hubble constant Hq — 62.3 ± 
1.3(statistical) ± 5.0(systematic). So considering the Sandage calibration, we have to 
add 0.02mag to distance moduli of the Gold sample and substract O.Olmag from the 
values given in the ESSENGE sets. 
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Table 1: Number of SNe indicating acceleration or deceleration at 95% 
and 99% CL for the different data sets and calibrations. Also given is the 
total number of SNe in each set. The most different results are highlighted. 





Gold 2007 (MLCS2k2) ESSENCE (MLCS2k2) ESSENCE (SALT) 




Riess 


Sandage 


Riess 


Sandage 


Riess 


Sandage 


acceleration (95% CL) 


37 


27 


37 


14 


96 


53 


acceleration (99% CL) 


13 


7 


11 


3 


64 


30 


deceleration (95% CL) 


2 


4 


3 


7 


1 


7 


deceleration (99% CL) 





1 


2 


3 





2 


number of SNe 


182 


182 


162 


162 


178 


178 



4. Results for a fiat universe 

4.1. Single SNe 

The relevant parameter for our analysis is A/Xj given by equation (7). Its values are 
plotted in figure 1 for the SNe of the three data sets and the two different calibration 
methods. The curve for a flat ACDM model (with Q,^ — 0.3) is shown in the plots 
in order to give the reader a notion of how the redshift-dependency of A// could look 
like. Most of the data values are positive which indicates an accelerated expansion. A 
similar diagram has been presented by Gong et al. [17] for single SNe as well as for 
binned SN data where they used a combined Gold and ESSENCE data set. They kept 
the arbitrary value of the absolute magnitude M given in the set and then determined 
the Hubble constant by using SNe with z < 0.1 to he Hq — 66.04. At this point they 
stopped their analysis by concluding that due to the large number of SNe that he above 
the curve of a universe with q — 0, accelerated expansion is evident. But the fact that 
there are also several data points below that curve and that the errors (Tj (with a typical 
value of 0.23mag for MLCS2k2 and O.lSmag for SALT) are approximately of the same 
size as the values Afii give rise to the question how certain acceleration really is. Thus, 
we will provide a more quantitative analysis in this work. 

We counted the SNe of each set whose values of A/Ji were above the action limit 
y4a(95%) = 1.645(jj and those above Aa(99%) = 2.326ai (which indicates acceleration at 
a 95% and a 99% CL, respectively) and those with values below Ad{95%) — — 1.645(7i 
and below Ad(99%) = — 2.326(7^ (which indicates deceleration). The results are shown 
in table 1. The clearest evidence for an accelerated expansion is given by the ESSENCE 
(SALT) set with the Riess cahbration with 64 SNe indicating acceleration at a 99% 
confidence level and none indicating deceleration. Also in the other sets accelerated 
expansion is preferred with the exception of ESSENCE (MLCS2k2) in the Sandage 
calibration where none of the two expansion histories is preferred in the test using the 
action limit A{99%). 

It is noticeable that the two light-curve fitting methods for the ESSENCE set 
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(a) Gold 2007 (MLCS2k2), Riess calibration 



(b) Gold 2007 (MLCS2k2), Sandage calibration 
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(c) ESSENCE (MLCS2k2), Riess calibration (d) ESSENCE (MLCS2k2), Sandage calibration 
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(e) ESSENCE (SALT), Riess calibration 



(f) ESSENCE (SALT), Sandage calibration 



Figure 1: Magnitude A/Xj, as defined in equation (7), for the three data sets and the 
two cahbrations. SNe that indicate acceleration at 99% CL are plotted in red (x), those 
that indicate deceleration in blue (^). Also shown is the curve for a fiat ACDM model 
with Qrn = 0.3 (which is not a fit to the data). Note that the SALT method leads to a 
larger spread in A/Xj, whereas the Gold set extends to higher redshifts. 
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Figure 2: Differences of apparent magnitudes obtained by the SALT and the MLCS2k2 
fitting method (a) and these differences devided by (Xj (b). 



yield very different results. (For a discussion on these differences at low redshifts, see 
[22].) This could be an indication that at least one of the two fitting methods does 
not give the correct result. In order to consider this possibility we take a look at the 
differences between A/ij obtained by SALT and A/ij obtained by MLCS2k2. They are 
shown in figure 2a for the 153 SNe contained in both ESSENCE sets. For some SNe 
the difference is very large (namely up to one magnitude). In order to see if the two 
sets are consistent, we need to know how large the difference is in terms of the error 
ai = [af (SALT) + cT2(MLCS2k2)]^/^ The result is shown in figure 2b. We want the 
systematical error due to the different fitting methods to be smaller than the statistical 
error of the observational data. Thus, we discard all SNe with a difference in A/i, 
larger than lai and those that are only contained in one of the two ESSENCE sets, 
which leaves 129 SNe in the sets. Again we counted the SNe indicating acceleration or 
deceleration for each set and each calibration method. The result is shown in table 2. 
The outcome of this test did not change qualitatively by discarding suspicious SNe from 
the sets. Thus, we will use all the SNe of the ESSENCE sets in the following. 

The error of the distance modulus contains the peculiar velocities fpec of the SNe. 
In the ESSENCE sets ay is assumed to be 400km/s for all SNe. We wanted to know 
how sensitive our test is with respect to changes of the peculiar velocity. Table 3 shows 
that varying cr^ almost does not change the result with the exception of ESSENCE 
(MLCS2k2) in the Sandage calibration, where the number of SNe indicating deceleration 
at a 95% CL is 7 for ay = 400km/s and 500km/s, but 11 for ay = 300km/s. 

4-2. Averaging over SN data 

Until now, we have only made tests for single SNe. The problem in combining these 
data is that the quantity A/i depends on the redshift z. Thus the data from different 
redshifts do not have the same mean value. Nevertheless it is possible to average over 
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Table 2: Number of SNe indicating acceleration or deceleration for SNe 
of the ESSENCE sets with |A/ii(SALT) - A/ii(MLCS2k2)|/cTi < 1. The 
total number of SNe in each set is 129. Again, the most discrepant results 
are highhghted. 





ESSENCE (MLCS2k2) ESSENCE (SALT) 




Riess 


Sandage 


Riess 


Sandage 


acceleration (95% CL) 


33 


11 


69 


30 


acceleration (99% CL) 


8 


2 


40 


14 


deceleration (95% CL) 


1 


2 





6 


deceleration (99% CL) 





1 





1 



Table 3: Number of SNe indicating acceleration or deceleration for SNe 
of the ESSENCE sets with different pecuhar velocity dispersions (7„. 

ESSENCE (MLCS2k2) ESSENCE (SALT) 



(7y [km/s] 




SaiidMgc 


r!i(>ss 


Santiago 


300 acceleration (95% CL) 


37 


14 


99 


53 


acceleration (99% CL) 


11 


3 


64 


30 


deceleration (95% CL) 


3 


11 


1 


7 


deceleration (99% CL) 


2 


3 





2 


400 acceleration (95% CL) 


37 


14 


96 


53 


acceleration (99% CL) 


11 


3 


64 


30 


deceleration (95% CL) 


3 


7 


1 


7 


deceleration (99% CL) 


2 


3 





2 


500 acceleration (95% CL) 


36 


13 


95 


53 


acceleration (99% CL) 


10 


3 


63 


30 


deceleration (95% CL) 


3 


7 


1 


5 


deceleration (99% CL) 


2 


3 





2 


number of SNe 


162 


162 


178 


178 



A//j. Then the result of course depends on how many SNe from a certain redshift are 
used to calculate this mean value. Therefore, the average over all SNe of a set does not 
characterize the function Aii{z). But still, when combined with its standard deviation, 
the mean value can be an evidence for acceleration or deceleration. Thus, we define 

(10) 

where Qi — 1/crf. The mean value is calculated in such a way that data points with a 
small error are weighted more than those with large errors. The standard deviation of 
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(a) Riess calibration 



(b) Sandage calibration 



Figure 3: Magnitude A/x averaged over redshift bins of width 0.2 for different data sets. 



tfie mean value is calculated by 



a- 



(iV- 



(11) 



We start with averaging A/ij over redshift bins of width 0.2. For all data sets the 
averaged value of A/i increases with redshift (see figure 3) which could be expected for 
an accelerated expansion. The curve for a flat ACDM with = 0.3 has its maximum 
at z = 1.2 and then decreases again. Unfortunately, there are not enough SNe at 
high redshifts to state a possible decrease of A/j, at some point. The differences in the 
data points of the two ESSENCE sets again seem to be too large. MLCS2k2 gives 
smaller values than the SALT fitter for all redshift bins. As Afi{z) should be close to 
zero at low redshifts in a homogeneous and isotropic universe, the MLCS2k2 data look 
more consistent with our assumptions in the Riess calibration, whereas SALT gives the 
better results in the Sandage calibration. The Gold sample has sensible values in both 
calibrations. Afi devided by the error a-^ gives the evidence for acceleration in each 
redshift bin. As can be seen in table 4, the strongest evidence is given in redshifts 
between 0.4 and 0.8. 

Next we average over all SNe of each data set with a redshift z > 0.2. We discard 
SNe with a smaller redshift for the following reasons: (a) A/x(2;) is expected to be 
relatively close to zero for small redshifts. Therefore, nearby SNe do not contribute to 
the evidence for acceleration, (b) Considering only the nearby universe, local effects can 
modify the results for the distance modulus as the cosmological principal is not valid 
on small scales, (c) Another disadvantage of nearby SNe is that they were observed 
with many different telescopes and thus the obtained systematic error is potentially 
higher. Table 5 shows the mean values A/i and their standard deviations a-^ evaluated 
by using only SN data with z > 0.2. A/i is positive for all data sets. A/i divided by 
^A/I indicates the confidence level at which an accelerated expansion can be stated. The 
weakest evidence for acceleration is given for ESSENCE (MLCS2k2) in the Sandage 
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Table 4: Statistical evidence A/i/a-^ within the given redshift range for 
a flat and an open universe. 

Gold 2007 (MLCS2k2) ESSENCE (MLCS2k2) ESSENCE (SALT) 



z 


Riess 


Sandage 


Riess 


Sandage 


Riess 


Sandage 


flat universe 














0.0 - 0.2 


2.1 


0.4 


1.6 


-2.8 


5.7 


0.8 


0.2 - 0.4 


4.2 


2.9 


3.9 


0.7 


7.3 


4.3 


0.4 - 0.6 


9.8 


7.6 


7.1 


3.1 


12.1 


7.7 


0.6 - 0.8 


5.4 


4.0 


7.7 


3.7 


8.5 


4.8 


0.8 1.0 


4.8 


3.5 


6.8 


4.1 


6.0 


4.2 


1.0 - 1.2 


2.8 


2.2 










1.2 - 1.4 


2.3 


1.8 










open universe 














0.0 - 0.2 


2.1 


0.4 


1.5 


-2.8 


5.7 


0.7 


0.2 - 0.4 


3.5 


2.2 


3.2 


-0.0 


6.6 


3.6 


0.4 - 0.6 


7.4 


5.2 


5.5 


1.4 


10.2 


5.8 


0.6 - 0.8 


2.8 


1.4 


4.9 


0.9 


5.9 


2.2 


0.8 - 1.0 


1.4 


0.2 


3.9 


1.2 


4.0 


2.2 


1.0 - 1.2 


0.9 


0.4 










1.2 - 1.4 


-0.2 


-0.7 











Table 5: Mean values and standard deviations of A// obtained by using 
only SNe with z > 0.2 for a flat universe. 

Gold 2007 (MLCS2k2) ESSENCE (MLCS2k2) ESSENCE (SALT) 



z 

'Ajl 



Riess 


Sandage 


Riess 


Sandage 


Riess 


Sandage 


0.63 


0.63 


0.54 


0.54 


0.51 


0.51 


0.2196 


0.1655 


0.2398 


0.1056 


0.3457 


0.2115 


0.0167 


0.0167 


0.0201 


0.0201 


0.0203 


0.0203 


13.1 


9.9 


11.9 


5.2 


17.0 


10.4 



calibration. Here the mean value lies 5.2a above 0, i.e. above the value for a universe 
that neither accelerates nor decelerates. In the other cases the confidence level is even 
larger, up to 17. Oct for ESSENCE (SALT) in the Riess calibration. 

5. Open and closed universe 

Although there are good reasons to believe that the universe is flat we give up this 
assumption in the following section. In an open universe the luminosity distance of a 
universe that neither accelerates nor decelerates is given by 

dL(-2) = -^^sinh(v^ln(l + z)) (12) 
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Table 6: Statistical evidence A/i/a-^ for an open universe (obtained by 
using SNe within the redshift range 0.2 < z < 1.2), a flat and a closed 
universe (0.2 < z). 

Gold 2007 (MLCS2k2) ESSENCE (MLCS2k2) ESSENCE (SALT) 
Riess Sandage Riess Sandage Riess Sandage 



open universe 


8.0 


4.9 


8.8 


1.8 


13.8 


7.2 


flat universe 


13.1 


9.9 


11.9 


5.2 


17.0 


10.4 


closed universe 


13.1 


9.9 


11.9 


5.2 


17.0 


10.4 



and thus depends on the density parameter of the scalar curvature fife, increases 
with increasing fife and thus the evidence for acceleration becomes weaker. As we are 
interested in the lower limit of this evidence we have to take the highest possible value 
for the scalar curvature, i.e. fi^ = 1, which corresponds to an empty universe. (Here we 
allow ourself to make use of the Einstein equation.) Then equation (7) for A/ij changes 
to 

AiJ,i ^ Hi - i_i{q ^ 0) ^ Hi - 5 log 



,1 + Zi) sinh[ln(l + Zi)] 

-no 



-25. (13) 



The evidence for accelerated expansion is then calculated in the same way as for a 
fiat universe. The result obtained by averaging over redshift bins is shown in table 4. 
For the overall average we only used SNe between redshift 0.2 and 1.2 because including 
higher redshifts would weaken the evidence. This is due to the fact that the values 
of A//j become negative when the phase of deceleration within the redshift range over 
which is integrated is large enough. The result is shown in table 6. As could be expected, 
the evidence is now much weaker than for a flat universe. But we still find a hint of 
acceleration at I.80". 

For a closed universe we have a different situation: Here c/l decreases with increasing 
spatial curvature. Thus the lower limit of the evidence for acceleration is given in the 
case of the lowest possible curvature which corresponds to a flat universe. 



6. Conclusion 



We tested the cosmic expansion without specifying any density parameters or 
parameterizing kinematical quantities. This was not possible without assuming a certain 
calibration of M and Hq. Therefore, we considered two very different calibrations. We 
also considered two different data sets and two different light curve fitters and varied 
the peculiar velocity dispersion. Using single SNe for the test has already given a 
clear indication of acceleration in a flat universe although a few SNe strongly favour 
deceleration. We find large systematic effects already at the level of individual SNe. 

Similar analyses have already been done by other groups [16, 17]. But note that the 
work of Santos et al. [16] exhibits two major shortcomings. First, they did not cahbrate 
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the SN data consistently. They took a certain value for Hq but kept the arbitrary 
value of the absolute magnitude M given in the set which led to distance moduli that 
are too high. Thus, they concluded erroneously that there was a recent phase of super- 
acceleration. The second problem is that they did not realize that due to the integration 
over redshift this method is not suitable to determine the transition redshift Zt between 
deceleration and acceleration. Thus, their value of Zt does not hold. The shortcoming 
of the work of Gong et al. [17] is that they have not applied any statistics at all. 

Instead of only considering single SNe, we obtained a more significant result by 
averaging over the SN data of each set. Here it is justified to discard nearby SNe in order 
to decrease systematics from using different telescopes or from local effects. We argue 
that we reject the null hypothesis of no acceleration for a spatially flat, homogeneous and 
isotropic universe at high confidence (> 5a). However, the results show large differences 
for each data set, calibration and light-curve fitter. E.g. changing the calibration from 
Sandage to Riess calibration in a flat universe using ESSENCE(MLCS2k2) increases the 
evidence for acceleration from 5.2a to ll.Qcr. For the two different light-curve fitters 
applied to the ESSENCE SNe using the Riess cahbration, we get values of ll.Qcr and 
17.0(7. 

Wood-Vasey et al. argue that the differences due to different fight curve fitters are 
not important as they disappear when marginalization is applied [6]. This is true if the 
data are only used to determine the parameter values of a certain cosmological model 
as in this case no calibration of M and Hq is needed. But nevertheless, the fitters 
should give the same result for a given absolute magnitude M. If this is not the case, a 
systematic error in the determination of the apparent magnitude m is introduced by at 
least one of the fitters. Systematics due to different fitters have also been described in 
[22] and those within the Gold sample in [23, 24, 25]. 

The different results obtained by the two calibrations are very surprising. Our 
test only depends on the reduced distance modulus M = M — 5 log Hq + 25 and not 
on the absolute magnitude or the Hubble constant individually. (Note that this is 
the quantity which is marginalized when fitting cosmological parameters.) Starting 
from different calibrations Mi and M2 leads to different values of Hq, determined by 
observation of SNe, in such a way that Ai should in principle be the same for both 
values of M. Thus our test should not depend on the calibration. However, the fact 
that the calibration changes our result significantly can be explained as follows: Hq 
is determined by observing SNe. Riess et al. and Sandage et al. use different sets of 
SNe and different fitting methods for this determination. The systematic errors due to 
different sets and fitters then of course also infiuence the result for the Hubble constant 
leading to different values M.i and M.2- 

For a conservative conclusion, we need to take the set that gives the weakest 
evidence for acceleration, namely ESSENCE (MLCS2k2) in the Sandage calibration. 
Using this set, accelerated expansion can be stated at 5a if we assume a spatially flat or 
closed universe. In an open universe the evidence is much weaker, namely 1.8a. These 
results only hold if the correct analysis of SNe is somewhere in the range of the cases we 



How strong is the evidence for accelerated expansion? 



14 



considered here. But as we observe enormous systematic effects in our test, we cannot 
be sure of tliis assumption. Tlius, it is a major issue to better understand liow SNe fiave 
to be observed and analyzed. 

Remember tliat results of this work are only valid if the universe is homogeneous and 
isotropic. The situation changes dramatically if we give up those assumptions. There 
exist enormously large structures in the universe, an example is the Sloan Great Wall 
with an extension of ~400 Mpc [26]. Besides, there are also big voids and superclusters 
on a 100 Mpc scale. This roughness of the universe could spoil the validity of equation 
(1). Indeed it is possible to construct inhomogeneous models that can describe the 
observational data without the need of an average acceleration [27, 28, 29, 30]. 

Some observational evidence for inhomogeneity and anisotropy in SN Hubble 
diagrams has recently been presented by [31, 32], probably due to large scale bulk 
motion and perhaps systematic effects. A non-trivial scale dependence of the Hubble 
rate due to the so-called cosmological backreaction [33, 34, 35, 36, 37, 38] has been 
shown in [39]. This effect of averaging can also mimic curvature effects [40] and thus 
strongly influence the reconstruction of the dark energy equation of state, as has been 
shown in [41]. 

Thus we think a major task must be to establish the acceleration of the universe 
independently of the assumption of strict homogeneity and isotropy. 
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